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A logical foundation of equilibrium state density functional theory in a Kohn-Sham type formula- 
tion is presented on the basis of Mermin's treatment of the grand canonical state, it is simpler and 
more satisfactory compared to the usual derivation of ground state theory, and free of remaining 
open points of the latter. It may in particular be relevant with respect to cases of spontaneous 
symmetry breaking like non-coUinear magnetism and orbital order. 
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I. INTRODUCTION 

Modern ground state density functional theory (DFT) 
for an inhomogeneous system of identical particles, hav- 
ing early roots in the work of Thomas and Fermi, was 
pioneered by the seminal papers by Hohenberg, Kohn 
and Sham.— 1^ It was later generalized by the constrained 
search concept of Levy"^ and finally put on a mathemat- 
ically rigorous basis of functional Legendre transforms 
by Lieb.- Meanwhile, DFT for the quasi-particle self- 
energy labeled by acronyms like GW or LDA-I-DMFT, 
and time-dependent DFT for dynamical processes, re- 
lated to Keldysh Green's functions, where developed, and 
all these theories layed the ground for enormously suc- 
cessful model approaches (by use of model functionals) 
to simulation of molecules and solids of any complexity. 

A generalization of the ground state DFT to thermo- 
dynamic equilibrium states by Mermin^ which appeared 
shortly after the work of Hohenberg, Kohn and Sham, 
has been mentioned from time to time^^ but otherwise 
has been largely ignored so far. 

As it sometimes happens^ a formalism for tempera- 
ture T = need not be equivalent to that for T J, 
which latter case is always the relevant case in physics. 
Although ground state DFT is largely settled now, some 
uneasy feeling remains in connection with the density 
functional failing to be differentiable in some cases (where 
nss' I— > Vss' is not unique), notably in spin DFT.^'^° In 
Ref. 1^ it is argued that these problems may reduce to the 
ordinary well understood gap problem, now for the spin 
subsystems separately, if one restricts consideration to 
homogeneous external magnetic fields only. This might 
seem a reasonable restriction since applied fields in labo- 
ratory can hardly vary over microscopic distances. How- 
ever, in the very topical cases of spontaneous symmetry 
breaking with respect to the interplay of non-coUinear 
magnetism with orbital order, in a statistical treatment 
one has to resort to the trick of Bogolubov's quasi-means 
by applying a suitable infinitesimal symmetry-breaking 
external field, otherwise statistical ensembles would not 
reproduce the broken symmetry. In the just mentioned 
cases this implies a microscopically inhomogeneous field, 
and one would like to rely on a situation where every- 
thing is fine at least in an infinitesimal vicinity of such 
a symmetry-breaking field. The good news is that this 



is indeed the case and the needed functional derivatives 
always exist for T > 0. This is shown with the help of 
Mermin's approach in the sequel. DFT is a rigorous the- 
ory for volume < oo and for temperature T > 0. For 
V = +00 a ground state wave function may not exist 
and for T = the functional derivative of the density 
functional may not exist. The theory then may be ap- 
plied for V ^ oo (as is routinely done with refinement of 
the grid in fc-space; a discrete regular k grid means pe- 
riodic boundary conditions with a finite periodicity vol- 
ume) and for T 4- 0. (Contrary to the case of adiabatic 
molecular dynamics^ the temperature of a finite system 
in an equilibrium state has a well defined meaning in the 
average over states the system may be in after in had been 
for a long time in contact with a large thermal bath.) 



II. SUBTLETIES OF 
HOHENBERG-KOHN-SHAM THEORY 

Originallyl DFT was built for systems in an (arbitrar- 
ily large) box of finite volume which conveniently can 
be replaced by periodic boundary conditions meaning to 
treat the position space of the particles as a three-torus 
T"^ of finite volume for instance x = x + L, y = 

y + L, z = z + L, |T3| = 

Let the Hamiltonian be 



H = T + W + V, 



^2 

f = — / v4>'i{x)V4>{x)dx, 



(1) 



(2) 



W= I i;\x)i^\x')w{r,r')i){x')i){x)dx'dx, (3) 



V— 'ii)'^{r,s)v{r)'ij^{r,s)dx, 



(4) 



where '^l^ix), x — (r, s) with particle position r and spin 
variable s, is the field operator of the particle field and 
Jdx^EsJd'r. 

Then, given any particle number N, a normalized 
ground state (GS) wavefunction (WF) ^(a;i, . . . , xn) and 
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a GS density n{r) exist for any reasonable external po- 
tential v{r) and for any non- negative pair interaction 
w{ri,rj). As is now standard,— one allows for all poten- 
tials with the only condition that Jj3 \v\'^^'^ cPr < oo, that 
is, V G L^/^(T^). Potentials of arrays of finitely many 
point charges in the belong to this space^^ii, and the 
Hamiltonian — T + V oi interaction-free fermions is 
bounded below for any such potential. Then, this also 
holds true for Hamiltonians ([Ij, if w{r,r') > 0. Since 
the space has finite volume, all considered Hamil- 
tonians have discrete spectra with at most finite de- 
grees of level degeneracy. (Lieb allowed the position 
space to be the real vector space R'^ of infinite volume 
which caused many problems with the continuous part 
of the spectrum of Hamiltonians. He then had to restrict 
n G L'^(M^)nX^(R^) since the density must integrate to a 
finite particle number TV over the infinite space R'^. This 
led him allow for potentials v e L^/^{R^) + L°°{R^). In 
the three-torus every function n G L'^{T'^) may be nor- 
malized to integrate to a given N.) 

The lemma by Hohenberg and Kohn^ states the unique 
mapping n H> w from ground state densities (degenerate 
ground states allowed'*) to external potentials, on which 
basis the Hohenberg-Kohn density functional 



-fHK["-] — E[v[n\,N] — {n,v[n\) 



(5) 



is defined for any ground state density n G An- 

An = {n comming from an iV-particle GS-WF}. (6) 

In ^ v[n] means the potential causing a WF-GS density 
n, and henceforth we use the notation of linear function- 
al 



{n,v) 



T3 



n{r)v{r) d^r. 



(7) 



As is easily seen, {n, v) cancels an equal term in the GS 
energy E[v[n], N], so that Fhk does not any more depend 
on v[n\. The functional F = Fhk might be used in the 
variational principle by Hohenberg and Kohn 

E[v, N] = min{i^[n] + {v, n)\{l,n) ^ N} (8) 

n 

where means a set of elements A with property 

B. The crucial point for the possibility to solve this 
problem with the help of Euler's equation is the know- 
ledge of the variational domain for n and the existence 
of the functional derivative of F. Would the derivative 
of F at the minimizing density nmin exist, it would be 
5F[ni^{i^/ 5n = —v + /i, where /i is the Lagrange multi- 
plier for the constraint in ([8]). For i^HK, unfortunately 
neither the domain of definition An '^^ explicitly known 
nor is anything known about the existence of the func- 
tional derivative. We know that An C L'^iT'^), but Lieb 
has shown^ that An is not convex: There are densities 
n ~ '^j^CiUi, Ci > 0, J^i'^i — 1' which are not in An 
while the Ui all are in ^^v- This is, why nowadays more 
general definitions of F[n] are used. 



As the theory can be build for any reasonable pair in- 
teraction w, the interaction- free case u; = is of some 
help. Further on this case will be denoted by a super- 
script 0, but the corresponding density functionals F^ 
will as usually be denoted by T since they obviously re- 
duce to the kinetic energy of an interaction-free system 
with GS density n. In this case, an alternative to ^ is 
the density matrix (DM) functional 



mm 

{Pk,Vk} 
0<pk<l 



^Pk{(pk\i\ipk) 



k 



(9) 



with i = —h^/\/2m. For any A^-particle DM state (not 
only GS), 



T[{pk,'~Pk}] = ^Pk{^Pk\i\ 



Vk) 



<Pk < 1, ^Pk = N, {ipk\'Pk') = 5kk' 



(10) 



is the general expression of the kinetic energy (with the 
set of orthonormal orbitals depending on the state). 
Given a potential u, the orbitals with {t -t- v)Lpk = (fk^k 
and occupation numbers pk = 1 for Sk < £Nt Pk = for 
Sk > £n minimize (O for the corresponding GS density 
n G Ay)j^n- The GS is unique if the highest occupied 
level Eat is not degenerate. The GSs and their energies are 
obtained from the Kohn-Sham (KS) variational principle 



min r[{pfc,V3fc}] V((/3fc|-y|v?fc) ^ 

Pk.V>k} i , , ) 

^1 ^ ^, — AT ft—L 



0<Pk<l,J2pk=N 
{Vk\Vk')=^kk' 

(11) 

while ^ is defined for the density n of any A^-particle 
state, that is,^ on the domain 

Jj^^{n \ n{r) > 0, V{n^^^) G ^^(T^), (1, n) = N} 

(12) 

which is a convex subset of L^{T^). (It has been showni 
that a minimum (|9]) exists for every n G Jn-) 

Let for the sake of simplicity v have a non-degenerate 
GS (which is a single determinant of orbitals in this case) 
with density n. It is easily seen that the most gen- 
eral variation permitted by ([9]) is a linear combination of 
{Sipk = Xk(pi^}, Sk < £jv, Eife > ejv, Afc 0. To lowest 
order in the Xk it yields 5Tdm = 2Re J2k ^ki^h \'t\^k) = 
2Re '^k^k(}Pik\'v\^k) = {v,Sn). These variations reach 
every single determinant state in a neighborhood of 
the considered GS (with respect to the iif^(T'^^)-norm 
||^i|j2 ^ /(|*|2 + |V^'p)d^^r). Since^ single determi- 
nant states map onto Jn , the corresponding Sn is a gen- 
eral variation in a neighborhood of n in J7jv (relative to 
the L'^(T'^)-norm), and hence the functional derivative 
of the convex functional Tdm exists at n as a (Frechet) 
derivative in Jn and equals v. The argument can be 
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generalized to the case of a degenerate GS, that is, for 
all n € -^DM.TV' 

Now, let n £ Jjy \ ^. Such densities exist, for 
instance densities having nodes cannot be in A^j^ jy- 
sume that the derivative of Tdm exists at that n. This 
means that there is some u with STdm = (m, 6n) for all 
permitted 6n. Since Tdm was shown^ to be convex, the 
assumption implies that n minimizes TDM[fi'] — {u,n') 
and hence is a GS density to the potential —u in contra- 
diction to the presupposition. Tdm has nowhere outside 
of A!^yi a functional derivative. 

KS theory in the interacting case w > now uses the 
splitting 



F[n] = TdmN + En[n] + E^cH. 



(13) 



which defines the density functional Exc[n] through the 
preceding ones. While this definition is correct, nothing 
can be said on the existence of the functional derivative of 
-E'xcN for GS densities n S .4dm, Af since we do not know 
the sets A^y^ n ^^'^ A'dm,n and cannot assume .Adm.tv C 
-^DM Af- Even though F has a functional derivative for 
n G .4dm, Af (see below), -Exc["-] can only have one there 
if Tdm has one. 

Like in with n — X^fePfelVfcP a KS variational 

principle is set out with the KS equation as the corre- 
sponding Euler equation. Now one may assume 



TdmM-I-SxcN 



0<pfc<l 



N 



min <^ V'pfc(v5fe|iVfc) + ^^xc[{Pfe,¥'fc}] 

{Pk,Vk} ' 



k=l 



(14) 



AT 



fc=i 



leading to a KS equation with a (non-linear) exchange 
and correlation potential operator 



1 5£' 



xc 



Pk Siplir) 



vxcVkir), (vxcVk)* = vxcV't , (15) 



if one assumes that Sxc depends on pk and on a Hermi- 
tian form of the (pk only. Like in ^ , the right hand side 
of (|14p might exist (not proven so far; only £x[{Pk, ^k}] is 
a simple explicitly known expression). By inserting (jl4p 
together with n — ^f.P'^l'Pkl'^ C3) it is seen that 
£xc (if it exists) has derivatives with respect to pk and 
ipk since all the other terms in the equation have them 
(even where the left hand side of p4|) has no derivative 
with respect to n). Hence, -Cxc will in general not be 
a local potential function, it may in particular be or- 
bital dependent (vx is non-local, its orbital dependence 
is canceled by inclusion of the orbital dependent self- 
interaction in both vn and vx] the model XC poten- 
tials with partial self-interaction correction or in LDA+U 
models are non-local and orbital dependent). Would vxc 
exist as a local potential, then the KS equation would al- 
ways yield a solution n S A'^^ ^ and hence there would 



be .4dm, Af C .4dm n which can by no means be taken for 
granted. 

The only density functional F for which the issue of the 
existence of the functional derivative can be addressed in 
general is the Legcndrc transform'* 

F[n]= sup {E[v,N]- {n,v)} (16) 

vex* 

for both cases, w = and w ^ 0. It is convex and de- 
fined on the whole functional space X (it may take on 
the value -|-c» in part of X), and if, given n, there exists 
a unique maximizing v, then this is the functional deriva- 
tive of F[n]. Since v is indeed up to a constant uniquely 
determined by any GS n, the functional derivative of this 
F exists at least for n G .4dm, at as a derivative ('gradi- 
ent', more precisely Frechet derivative) in the hyperplane 
{neX\{l,n) =iV}. 

Less clear is the situation in spin DFT.^ Now, also F 
need not have a derivative for GS densities. 



III. A FEW ESSENTIALS ON LEGENDRE 
TRANSFORMS 

Let X — X** and X* be two mutually dual functional 
spaces, that is, X* comprises all norm-continuous lin- 
ear functions on X and vice versa. ((K.^)* = is 
the space of all gradient vectors to functions on M^; 
(i3(T3))* = i3/2(T3) and vice versa.) The Legendre 
transform f*{n), n E X of a, function f{v), v € X* is 
defined as 



f*{n) ^ sup {{n,v) - f{v)}. 
vex* 

A second Legendre transformation yields 

f**{v)^sup{{v,n)~-f*in)}. 
nex 



(17) 



(18) 



All we need is 

1. f*{n) is a convex function of n, no matter what 
f{v) is; if f{v) is convex, then f**{v) — f{v); in 
general f**{v) < f{v). 

2- f{v) + f*(n) > {v,n); if, for convex /, f{v) + 
f*(n) = {v,n), then v G df*{n) and n G df{v). 

In the second statement df{v) is the subdifferential of 
the convex function / at point v: the set of all linear 
functions {n,v') so that f{v') > f{v) + (n, {v' — v)) for 
all v' G X* . If this set consists of a single linear function 
only, then this linear function is the (total) differential 
df{v), that is, n is the derivative of / at v. 

To elucidate these properties one may consider con- 
vex functions of one real variable, f{N) and /*(//) (see 
e.g. Fig. 11 of Ref. [111). Put a supporting tangent to 
the graph of / at point N (a line having the common 
point {N, f{N)) with the graph of the function and be- 
ing nowhere above). The tangent has a slope fi. The 
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sign carrying distance from the intersection point of this 
hne with the /-axis to the coordinate origin is If 
/ has a derivative at N, then its value is /i. It is easily 
seen that, if the derivative of / jumps at N, then there 
is a (closed) interval [/ii,/i2] from the left derivative 
to the right derivative /i2 i^J■l may be —oo or fi2 may be 
+00), and f*{^i) is linear on this interval, the interval 
being the subdifferential df{N). Inversely, if the convex 
function / is not strictly convex, but has a linear depen- 
dence on some interval with slope ^, then the derivative 
of /* jumps at that /i. 

This simple geometric picture readily transfers to the 
general case: take a tangent hyperplane f{no) + {v, (n — 
no)) supporting the graph of f{n) at some no. The dis- 
tance from its intersection point with the /-axis to the 
origin is /*(«)■ If the derivative of / jumps at some n 
(and hence does not exist there), then there is a convex 
domain in w-space on which f*{v) is linearly depending 
on V, and vice versa. 

If the GS wave function is independent of some poten- 
tial change Sv called a 'phantom' potential perturbation 
in Ref. lOi, then the GS density n does also not change and 
the GS energy has a linear dependence const. + (Sv,n). 
Consequently, the functional derivative of F[n\ defined 
by (fT6|) does not exist at that n. This is precisely the 
role of 'phantom' potential perturbations in DFT. 



In the following tr will always mean the trace in the Fock 
space of the ip- Also, the natural abbreviation 

J2 J (vss'ir) ~ fiSss')ns'sir)d^r = {{v - n),n) (22) 



will be used. 

Now, fix the particle interaction w and, following 
Mermin^ (we try carefully to trace functional depen- 
dences and in doing so slightly deviate from Mermin's 
notation), consider for various external potentials v the 
functionals 



1 



ftyip] ^tr p[ H - pN + - Inp 



(23) 



As easily seen by direct substitution of (pn|. the grand 
canonical potential flp[v — p] is obtained as 

= -ilntre-^(^-'^^) =17,[p^[i;-m]]. (24) 

Moreover, as shown in Ref. [^, for any p > 0, trp = 1, it 
holds that 

n^ip] > n^ippiv-p]] ^ nfs[v-p] for p ^ ppiv-p]. (25) 



IV. UNIQUE MAPPINGS FOR T > 

We now move to temperature T > and to grand 
canonical states. We also generalize to spin DFT and 
allow for external magnetic fields coupling to the parti- 
cle spin but not to its charge (diamagnetic couplings as 
usually in spin DFT are neglected). Consider a system 
of identical particles in an external field Vss'ii")- Let the 
system be confined in a large box, or, placed in a large 
three-torus equivalent to periodic boundary conditions 
(regular fc-grid). Let the Hamiltonian be that of ([T][3]), 
but (HI) generalized to 



^ = i'Hr,s)vss'{rMr,s')d^ 



(19) 



The particle number operator is TV = J ip''^ {x)4'{x) dx 
so that H — pN depends on the combination w — /i = 
Vss' (r) - pSss> only. 

Fix the temperature /? ~ \/kT^ the chemical potential 
p and the external potential v. Then, the grand canonical 
state is 



p/3b - lA^ 



(20) 



If p > 0, trp = 1, is any state (density matrix), then 
iv p{V — pN) = J{v — p)n[p\ dx with the particle (spin) 
density 



In Mermin's approach, this inequality replaces the corre- 
sponding ground state energy property. It immediately 
implies that Q,ij[v — p] = miTip^y[p] is concave in v by 
the simple reasoning (we write in short Vi for Vi — Pi) 

Vlplavi + {1 — a)v2\ = 



min tr p aHi + (1 — a)H2 H In p > 

P \ P 



> a min tr pi Hi H In pi 

Pi V P 



n[p] = Uss'ir) = tT ptp(r,s)tp^{r,s'). (21) 



(26) 



(1 - a) min tr P2 + 7^ In P2 = 
P2 \ P J 

= afifiivi] + {I - a)np[v2], 0<a<l, 

because a joint minimum of a sum cannot be below the 
sum of the independent minima of the items. 

As another advantage over the standard zero temper- 
ature theory, it follows immediately from ([50]) and (PT|) 
that the mappings {v — p) ^ pp ^ n are unique. There 
is no problem with degenerate states since degenerate 
states automatically get equal statistical weight in pp of 
(1201) . However, as usual spontaneous symmetry breaking 
is not covered by this statistical approach; it has to be 
treated by an infinitesimal symmetry breaking external 
potential v in the spirit of Bogolubov's quasi-means in 
Statistical Physics. Nevertheless, by virtue of (P5|) which 
also holds in the spin case, in the standard way Mer- 
min proved the analogue of the Hohenberg-Kohn lemma: 
n ^ (v — p) is unique for any n coming from a grand 
canonical ensemble at temperature l/k/3. In summary. 
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there are the unique mappings 

{v — fj.) < n 



grand canonical pp 



(27) 



On the functional domain (which may depend on W) 
Dfj — {n coming from some pp, /3 fixed} (28) 

one can write (v — p) p[n\ and also Hjsln] and pp[n\, as well 
as np[v — p] = n[pi3[v — p]] on the domain of admissible 
potentials v. (Denoting the distinct functions pf;[n\ and 
pplv] by the same symbol pp will cause no confusion.) 

Moreover, from the unique dependence of v — p on 
Pp it follows now also that VL^, [p\ for different Vi — pi 
is minimized by different pi, and hence the inequality in 
(|26p is sharpened into a strict inequality: ri^[w] is strictly 
concave. If equality would hold in (|26p. this would imply 
that the minimizing p is also a minimizing pi for vi and a 
minimizing p2 for V2 ■ This is the principal difference from 
the T = theory where E[v,N] is not always strictly 
concave in w^s' and is never strictly convex in N. 

V. THE DENSITY FUNCTIONAL 

As was already said, for electron systems it is well jus- 
tified to allow for all potentials 



(29) 



for which the integral J^^ \v — p\'^^'^ (Pr over the three- 
torus (of finite volume) is finite. Recall that the Hamilto- 
nian = T + V of interaction- free fermions is bounded 
below for any such potential, that this also holds true for 
Hamiltonians p1l51 [TO)) , if w{r,r') > 0, and that, since 
the space T'^ has finite volume, all considered Hamiltoni- 
ans have discrete spectra with at most finite degrees of 
level degeneracy. Then, f2^[w — p] of (IM)) is well defined 
on X* and smooth in the norm topology. 

In view of the concavity of ri^[t'], introduce the Leg- 
endre transform^'-'^^ of as Fp[—n]: 



Fp[n] = Fp[-n] = sup|-(n, t;) + ilp[v]^ 



(30) 



which as a Legendre transform is a convex functional of 
—n (or likewise of n), the dual variable to v. n G X** — 
X = L^{T^). 

Since the functional space X is reflexive, ^^(T^) = 
{L^{T^))**, the Legendre back transformation from (PO)) . 
-flplv] = sup„{-(u,n) - Fpi-n]} or, 

np[v]=ini{F,3[7i] + {v,n)} (31) 

represents the generalized Hohenberg-Kohn theorem 
(where equality holds since ri/3[w] is concave in v). The 



chemical potential p is further on put to zero which sim- 
ply means that single particle energies and potentials are 
measured from the chemical potential. 

For any density n £ Dp from (|28p . in analogy to the 
original Hohenberg-Kohn functional one may define 



[n] = [pp [n]]- {v,n), n e Dp 



and. from (1^51). have 



(32) 



np[v] = min \^p[n] + {v,n)] (33) 

since in view of (j27p any n ^ n[v] refers to pp[n\ ^ 
pp[n[v]]. From (PT|) and ([55)1 one infers that Fp[n] — 
^p[n] for all n G Dp, and that the infimum of ([51]) 
is always a minimum with minimizing density 
*/3NH] + {v,'np[v]) = np[v] < Fp[np[v]] + {v,np[v]), 
hence <I'/3[n] < for n S Dp, and by interchanging 

the role of (|3ip and ([55]) in the argument the opposite 
inequality is obtained. 

Moreover, for — and Fp[—n] like in general for 
any pair of mutual Legendre transforms it holds that 
Fp[—n] — flp[v] = —{n,v) implies v € dFp[—n] and —n G 
d{—^lp[v]), where dFp[—n] means the subdifferential on 
Fp at point —n and d{~Qp[v]) means the subdifferential 

on —flp at point v. Since c?ri^[ti] = Intr (e~^'^c?y) as 
easily seen from the definition of f2/3[w], its first (Frechet) 
derivative (for finite volume IT-^j) always exists so that its 
subdifferential contains only this one 'gradient'. Now, the 
reasoning after ([55)1 yields Fa[n/3[w]]— rj^[D] = —{np[v],v) 
and hence 



snp 

Sv 



np [v\ 



(34) 



As np[v] is a one-one mapping X* ^ Dp, for n E Dp one 
has inversely vp[n] and Fp[n] — flp[vp[n]] — —{n,vp[n]) 
implying 



SFp 
Sn 



-vp [n\. 



n e D 



(35) 



Note that while (1551) holds for n e Dp C X, the deriva- 
tive d/Sn on the left hand side is taken in X, that is, for 
any Sn X with \\5n\\ small enough. 

From the strict concavity and continuous differentia- 
bility of il.p[v] the differentiability of Fp[n] at every point 
n G Dp follows, that is, at every density n thermodynam- 
ically corresponding to some v at temperature {kl3)~^. 
Like in the T = theory^, ^ yields that Fp[ n] jumps 
to -|-oo if n < for an x-domain of non-zero measure: 
Take w = c > for some domain where n < and 
V = 0(= p) everywhere else. This v is admissible for 
arbitrary large c and H[v] is bounded below for such a v. 
Hence, ^p[v] is also bounded below and, as easily seen, 
the supremum ([501) obtained for c — > oo to be -l-oo. 
Assume now that np{xo) = for some xq. Since any 
np G Dp is continuous in x (any solution of the many- 
particle Schrodinger equation is continuous), there is al- 
ways Sn G X, Sn{xo) > 0, so that np(x) — eSn{x) would 
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be negative in a neighborhood of a;o of non-zero measure 
for arbitrarily smah |e| and the functional derivative psp 
would not exist for that np. Thus, the result (|35|) also 
implies nij{x) > everywhere for T > 0. (See also next 
section.) 



VI. INTERACTION-FREE PARTICLES AND 
BEYOND 

As is well known from Statistical Physics^^, in a non- 
interacting particle system the particles in a single parti- 
cle quantum state \(pk) may be treated as an independent 
subsystem even of a quantum ensemble with exchange 
symmetry. The corresponding statistical fermionic state 
is 



Pk = |>(1 -Pk){\ + \Vk)Pk{Vk\ 



(36) 



with occupation pfc of the orbital ifk and |) as the vacuum 
state. Accordingly we define^ 



T0[pk,Vk] = tr pk (^+ ^Inpfc 

+ i (pk Inpfe + (1 - Pk) hi(l - Pk)) 
and the density functional 



(37) 



,PkWk\ 



Tp[n\ — min 

{Pk,Vk} 
0<pfc<l 
('Pk\Vk')=^kk' 

(38) 

where the minimum taken over all orthonormal orbitals 
and orbital occupations which yield a given n exists like 
in the GS case. 

Now, the grand canonical potential is 



17° [w] = min|T^[n] + {v,n) 



0<pfc<l 
{Vk\Vk')=^kk' 



^"^1 <^ y''7^)[Pfc,'^fc] + (w,n[pfe,(^fc]) 

Ok.Vk} , 



(39) 



where again the Lpk must be orthonormal. 

Variation of the iy9^ under the last constraint yields for 
the minimizing orbitals 93° 



and variation of the pk yields 

Pl{fi)=fp{el) ^ 



(40) 



(41) 



which is the correct result in this physically trivial case. 
For any v G X* the minimum of ((M)) does indeed exist, 
and the minimizing density is 

"^H = (r) = M4)vli^, ^)^Tir, s') (42) 



so that N = J2k fl^i^k) relates the average particle num- 
ber N to the chemical potential fj,. Only the occupation 
numbers depend on temperature {kj5)~^ and on the value 
of the chemical potential /i from which v and the e° are 
measured. 

Here, n° > everywhere is intuitively clear because 
fp{el) > for all k. 

For densities minimizing p9p it obviously holds that 



n € D% 



(43) 



{Tfj[n\ replaces the density matrix functional TDivii'T-] of 
the ground state theory, Eq. ©.) 

Accounting for the Coulomb interaction of the elec- 
trons in mean-field approximation simply means to re- 
place V in the above Schrodinger equation by v+v^ where 



v^{r) 



nps's'ir') 



dx' 



(44) 



contains self- interaction. Since** ^^Pk G i^(T^) implies 
|(^fcp £ X^(T'^), taken as a KS ansatz n — X^fePfel'^feP — 
is sufficiently general for the density of an interacting 
system too. Densities of this type apparently form a 
convex domain D of the functional space X = L'^iT^) on 
which Tp[n] is also defined by (pS)) . [Jn C D for every 
real TV, < iV < cx).) 

AsDpdD also for ^ 0, by 



F0[n]=Tp[n] 



nss{r)ns's'{r') 



dxdx 



(45) 

for n e _D an exchange and correlation density functional 
Fg^^[n] is defined (since the other density functionals 
of this relation were previously defined or are explicitly 
given on D). Inserting here n = ^kPkWA"^ transforms 
(Plj) into a minimum search by varying Lp*^. and pk as 
above in the GS theory. The derivatives with respect to 
V?^ and Pk of Fp exist on the basis of ([35l) for {pk , } 
yielding n € D/^, and those of the second term on the 
right hand side of (HS)) are explicitly known. Hence, the 
situation with and F^'^ is like in the GS theory. We 
cannot expect Dp C D^. 

Formally, like in (|14p one may again assume 

Tp[n]+F^''[n] = 



mm 

{Pk,'Pk} , 
o<pfc<i 

{Vk\'Pk')=^kk' 



(46) 



which yields the KS equation 
with 



(47) 



1 f^F^C 

^eff^, = (v + v^ + v^^) ^k, v^'^'^k = - ^ , (48) 



Pk Sipl 
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and = fp from (|4T|) with £° replaced by Sk- Com- 

pare also the previous discussion of the property ([15]) of 
-xc_ Note that we did again not prove the existence of 
6F^'^ /Sn: v'^^, if it exists at all, need not exist as an 
orbital independent local potential, it might be non-local 
and orbital dependent. In this respect the situation is 
the same as for the ground state theory. 

Given the external potential v, the solutions of this KS 
equation determine, via the analogues of (|4T1 142|) without 
superscripts, the density n^lv] G minimizing the right 
hand side of (PT|) and hence providing the grand canonical 
potential 

riplv - ^] = Fplnplv - j {v - ^i)np[v - ^i\dx. (49) 

where we explicitly reinserted the chemical potential /i. 
The latter is related to the particle number N by 

- a^ = ^ = E//3(^'^-^) (50) 

k 



which is also confirmed by inserting (p5|) and the KS 
expression for n into (|49p . 

The whole theory, of course, as in the ground state 
variant again depends on the knowledge of the density 
functional [n\ and of T^'^ in the KS theory, both of 
which are hardly ever accessible (if the latter exists at all) 
and hence have to be modeled changing the exact theory 
into a model theory within a (nearly) rigorous frame. An 
early ad-hoc application is"'^. Since ((50|) rests on ([35]) . it 
can be used as a check for the quality of a model F^'^ [n] , 
for instance down to which temperature it can reasonably 
be used for a specific answer. 
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